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ABSTRACT 



Aims. To understand the dynamics of stellar interiors, we study the effect of rotation on turbulence. 

Methods. We consider turbulence induced by an arbitrary forcing and derive turbulence amplitude and turbulent transport coefficients 
(turbulent viscosity and diffusivity), first by using a quasi-linear theory and then by using a multi-scale renormalisation analysis. 
Results. With an isotropic forcing, the quasi-linear theory gives that the turbulent transport coefficients, both parallel and perpendicular to the 
rotation vector, have the asymptotic scaling fl -1 for rapid rotation (i.e. when the rotation rate CI is larger than the inverse of the correlation 
time of the forcing and the diffusion time), while the renormalisation analysis suggests a weaker dependence on Q., with fl~ 1/2 scaling. 
The turbulence amplitude is found to scale as Q° - Q. 1 in the rapid rotation limit depending on the property of the forcing. In the case of 
an anisotropic forcing with inhibited motion in the vertical direction, as should be relevant in a strongly stratified medium, we find that 
non-diffusive fluxes of angular momentum scale as Q. 1 - Q~ [ for rapid rotation, depending on the temporal correlation of the forcing. We 
discuss the implications of our result for the dynamics of stellar interiors. 
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1. Introduction 

Rotation and turbulence are ubiquitous features of many as- 
trophysical bodies and play a crucial role in the dynamics of 
these bodies. For instance, turbulence is very often invoked as a 
means to increase the tr ansport of angular mome ntum [e.g. the 
accretion rate in disks (Bal bus & Hawlevl 1 19981) 1 or the mix- 
ing of light elements in stars [e.g. to explain the surface de- 
pletion (see lPinsonneaultll 19971 and references therein)]. Since 
the rotation rate of stars varies with their mass (spectral types) 
and age, it is important to develop a general theory of turbulent 
transport which is valid for arbitrary rotation rate and other stel- 
lar parameters (such as molecular viscosity, particle diffusivity, 
etc). 

In the solar context, many authors have investigated the 
structure of turbulence in rotating bodies to explain the oc- 
currence of differential rotation in the convective zone. The 
main feature of this type of turbulence is the appearance of 
non-diffusive terms in the transport of angular momentum, 
which prevents a solid body rotation from being a solution o f 



the Reynolds equation ( Lebe dinskvll 1 94 1 ; Kippenhahnl l 19631) 



Starting from the Navier-Stokes equation, it is possible to show 
that these fluxes arise when there is a cause of anisotro py in the 
system , due to anisotropic background turbulence (see Riidiger 
19891 and references th erein) or inhomog eneities such as an 
underlying stratification dKichatinovlll987l) . To explain the ob- 
served differential rotation of the sun, it is also of prime impor- 
tance to understand the influence of rotation on the turbulent 



transport coefficients such as the turbulent diffusivity of par- 



ticles and the turbulent conductivity of heat. IKichatinov et al 



( 1994) have shown that for an isotropic original turbulence with 
long correlation time, all these coefficients scale as Q 1 (where 
O is the rotation rate of the sun) in the large rotation limit. 

In our recent publications, we have studied the dynam- 
ics of the Su n by takin g into account the prese nce of a 
strong shear dKirnl l2005t iLeprovost & Kim 2006) and the 
interaction of this sheared turbulence with different types 
of wav es that can be excited in the Sun due to magneti c 



fields (ILeprov ost & Kim 2007a; Kim & Leprovos t 2007a), 
stratification ( Kim & L eprovostl l2007bl) or global rotation 
(ILeprovost & Kimll2007bl) . The main purpose of these studies 



was to understand the dynamics of the thin radial shear laye r 
in the Sun, namely the solar tachocline ( Spiegel & Zahnl 1 992 ). 
where a sharp transition between latitudinal differential rota- 
tion in the convective envelope and nearly uniform rotation in 
the radiative interior takes place. 

In this paper, we consider the opposite limit where the 
differential rotation is negligible compared to the global ro- 
tation. This is relevant for the convection zone and also for 
rapidly rotating stars w hich are more massive than the Sun 
(Kii ker & Rudigen 120051) . Specifically, we investigate the ef- 
fect of rotation on the properties of turbulence; i.e. its intensity 
and its transport. First, using a quasi-linear theory of turbulence 
with an isotropic arbitrary forcing, we show that all the turbu- 
lent transport coefficients, both parallel and perpendicular to 
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the rotation vector, have the asymptotic scaling Q. 1 when the 
rotation rate is sufficiently large. The turbulence amplitude is 
however found to scale as Q° - Q -1 in the rapid rotation limit 
depending on the property of the forcing. We also consider the 
alternative case of an anisotropic turbulence with inhibited mo- 
tion in the vertical direction, as should be relevant, for exam- 
ple, in a strongly stratified medium. In this case, we find that 
the Reynolds stress does not vanish (unlike in the isotropic case 
where it vanishes) and that the transport of angular momentum 
scales as Q -2 — £2 in the rapid rotation limit. Secondly, we per- 
form a multi-scale renormalisation group analysis to improve 
on the quasi-linear theory. The result suggests that the trans- 
port scales as Q ' 2 , with a weaker dependence on f2 compared 
to that found by the quasi-linear theory. 

The rest of the paper is organised as follows: in Sect. [2] 
we solve the quasi-linear equations for the fluctuating velocity 
and density of particles in a rotating frame with an arbitrary 
external forcing. We then calculate the turbulent intensity and 
turbulent transport in the case of a homogeneous forcing with 
arbitrary temporal correlations in Sect. [3] In Sect. [Hand [gj we 
present results in the case of an isotropic and anisotropic forc- 
ing, respectively. Section|6]is devoted to the multi-scale renor- 
malisation group analysis. Finally, we summarise and discuss 
our results in Sect. [7] and provide the implications for the tur- 
bulence in stellar interiors in Sect. [8] 



2. Model (governing equations) 

Our starting point is the Navier-Stokes equation with a forcing 
term f combined with an advection-diffusion equation for the 
transport of chemical species, in a rotating frame with angular 
velocity fi: 



d,\ + v • Vv = -V> + vV 2 v + f - 2Q x v . 
V-v = 0, 
d,N + \-WN = DV 2 N , 



(1) 



where v is the viscosity of the fluid and D the diffusivity of 
particle. Note that, in the incompressible case considered here, 
the heat transport can be described by a n advection-diffusion 
equation (e.g. see IChandrasekharlll98ll) . Thus, our result for 
the turbulent diffusivity also holds for the turbulent conductiv- 
ity, provided that D is replaced by the thermal conductivity. 
In the following, we let = 2£l and assume (without loss 
of generality) that the rotation is around the x-axis with no 
large-scale velocity field (note that the com bined effect of ro- 
tation and large-scale shear is investigated in Leprovost & Kim 
2007bh . We then express the velocity field and the concentra- 
tion as the sum of a large-scale field and small-scale fluctua- 
tions: v = Up + u — u and N = No + n. In the quasi-linear 



theory dMoffattlll978l) . Eq. (Q]i can be linearised for the fluc- 
tuating fields and then Fourier-transformed to yield equations 
for the Fourier component X(k, t). To consider arbitrary Prandtl 
number P r = v/D, we introduce the following new variables 
X = exp[vk 2 t]X and X = exp[Dk 2 t]X to absorb the viscos- 
ity and diffusivity terms, respectively. Eliminating the pressure 
variable, the equations for the velocity field and concentration 



can be written: 
d 2 u x + (OqUx — 



d, 



hi(f) 



Oh 2 (f) 



(2) 



J + or 
j6a Oa „ h 2 (t) 

0,U, — O t U K -\ M V H , 

y y y 

u y — —(au x + f3u z ) , 
d,h = -u • ViVo , 

where a = k x /k y ,fi - kjk y , y = 1 +fi 2 and oj 2 - Cl 2 a 2 /(y + a 2 ); 
kit) = yfjt) ~ afy(t) ~Mz(f) and h 2 (t) = -f3f y (t) + f z (t). The 
homogeneous solution of the first equation can easily be found 
in terms of trigonometric functions. Using these solutions and 
the method of variation of constants, the solutions to the first 
three equations of Eq. (|2|i can then be derived: 



C' A hi?) . 0h 2 (t') 

u x = I dt { r-cos[wo(f —t)]+ sin[q)o(f - 

Jo \y + a z ^y + a 2 

u y = \ dt'i hl( '^ a (-cos[u) (t' -t)] 
Jo \y(y + a z )\ 

lit' - Ol) 



0] 



(3) 



0/3yjy + a 2 
H sin[(Wo( 



y^y + a 2 



^-6>sin[w (f' - t)] - 

Jo \yiy + a 2 )\ P 
Oyjy + a 2 



■ cos[w (f' - 



cos[wo(f' - f)] 



sin[w (f' - t)U ■ 



h 2 it')a 



rVr 



+ a z 



^6/3sm[ajo(t' - t)] - ^ 7 * ° cos[w (f' - f)]) 

where 6 is the sign of (£2a). 

In the following sections, we use Eq. ^ to calculate both 
the turbulet amplitude and transport (of chemicals or heat). To 
this end, we prescribe the correlation function of the forcing to 
be spatially homogeneous with a temporal correlation C(t): 

(kiku h)hj(k 2 , f 2 )> = (27r) 3 5(k! + fe) C(| t x - h\) 4> u (k 2 ) , (4) 

for i and j = 1 or 2. The functions <p t j are the power spectra of 
the forcing. 



3. Turbulent amplitude and transport 

3. 1 . Turbulent amplitude 

Using Eqs. ® and @J, it is straightforward to obtain the inten- 
sity of turbulence as: 



i(k) 



iy + a 2 ) 2 



(5) 



+A . C 2 0n(k) + AX 022(k) 



12 (r + a 2 F 2 22 (r + a 2 )/' 
2 ^ --(k) 

a 2 ) 2 

y 2a 2 <ft 12 (k) y a 2 (f> 22 ik) 



ilnf J \ n y 2 (y + c 



n y 2 iy + a 2 )V 2 22 y 2 (y + a 2 ) 
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cr. 



-f 

Jo 



e cos(wqt) C(t) dr . 



e - vk2r sin(w T)C(T) dr , 



(6) 



Here, the functions A'Jk, wo) are defined in Eq. jA.lj of These results in a general form are not shown here for the 
Appendix [A] they depend only on k and Q. and involve the sake of brevity and shall be presented in the case of isotropic 
temporal correlation function of the forcing C(t) only through and anisotropic forcing in Sect.|4]and|5] respectively, 
the functions k c and <x c , defined by: 

4. Isotropic forcing 

When the forcing f is isotropic, the correlation function of the 
forcing (with no helicity) can be written as: 

</KM)/;(k',0> = (2^) 3 C(|f-f'|) J F(fe)(5, 7 -^/P)5(k+k') .(12) 

The functions 0y in Eq. (01 are related to F(k) in ( fT2b by: 

0ii =y{y + a 2 )F{k) , 12 = and cf, 22 =yF(k). (13) 

Using Eq. ( fT3l ) in Eqs. (0, (© and ( [Tol l, we can obtain the tur- 
bulence amplitude and transport as follows: 

yF(k) k, 

■ ( // A' • 

(2: 



The turbulent amplitude in the z direction (u 2 ) can easily be 
obtained from (u 2 ), due to the symmetry of the problem, by 
making the following replacements: 



P 
hi 



faiP" 



afi~ 
hi ~ 



iP 



7 

-3 



yp 



-1 



(7) 



"22 



022 P 



-2 



3.2. Transport of angular momentum 

From Eqs. (f3]i and (01, we can also calculate the non-diagonal 
part of the correlation tensors (ujU j) to obtain the transport of 
angular momentum by turbulence: 

a0n(k) 

(2. 



In) 3 J 



(U X Uy) 



(2^) 3 J I 11 



r(r + a ) 



2>2 



(8) 



, 2fl 0i2(k) fl022(k) 



{UyU~) 



(2 



12 r(r + a 2 ) 3/2 22 7(7 + « 2 ) 

a 2 0n(k) 



2a 2 



7 2 (7 + a 2 ) 2 



2(k) 



12.. 2 



7 2 (7 + a 2 ) 3 / 2 



+ M: 



a 2 te(k) 



22 „2 



7 (7 + a ) 



Here, the functions Mfj and M z { . are given in Eq. ( IA.1I ) of ap- 
pendix[A]and depend only on A:, Qo and the functions k c and cr c 
defined in Eq. ©. 

3.3. Transport of particles 

Integrating the last equation of Eq. (01, we obtain the density 
of chemical species: 

hi = (-diNo) 



j Ui(f) dt' 
Jo 



(9) 



Note that the velocity to be integrated in Eq. (0 differs by an 
exponential factor from that given by Eq. (0. However it is 
more convenient to express the turbulent diffusivity of chemi- 
cal species Dj, defined by = -DjdjNo, in terms of u,. 
Specifically, we can obtain the following: 



(10) 



T (2n) 6 J 'J 

f dt' e (D - v)e S (ujikt, t')ui(k 2 , ?)> ■ 
Jo 

Using Eqs. (0-(0]i and following a calculation similar to that in 
Sect. 13. II we can calculate all the components of the turbulent 



diffusivity Dj. The results turn out to depend on the temporal 
correlation function C(r) through a new function: 



-f 

Jo 



e- DkT C{T)dT 



(11) 



<«v) 

(U X Uy) 
{UyU Z ) 
Dj 

DT 



d 3 k- 



(2tt) 3 J 7(7 + a 2 ) 



(7 + a 2 ) vk 2 
a 2 F(k) (y(J3 2 +a 2 ) 



(14) 



vk 2 a 2 



K e + 



0)Q 



v 2 k A + w 2 



(2n) 3 J 



3 aF(k) Kc 
dk , 

(7 + a 1 ) vk 1 

d\-^^Kc 



l - [d 3 k^ 

nf J (7 + 



(7 + a 2 ) vk 2 
yF(k) 



a 2 ) 



I(k, w ) , 



3 F(k)(J3 2 + a 2 ) 



I(k, w ) . 



7 (2tt) 3 J " " (7 + a 2 ) 

Here, the functions Xc, k c and cr e [defined in Eqs. (O and ( fTTT ll 
depend only on the modulus of the wave number; The function 
I(k, wo) is defined as: 



I(k,u) ) = 



1 



(£> - v) 2 /t 4 + w 2 
(D 2 - v 2 )fc 4 - w 2 



(15) 



(D + v) 2 k A + w, 



(«c + Q + 



2Dk 2 u) Q 



(D + y) 2 k A + w 2 

By performing the integration in the 7-direction and changing 
the colatitude variable to a> {) - |Qcos0|, we obtain the follow- 
ing result: 

Hni t / .,2^ 



(4) 



(u 2 ) = 



(2^) 2 |Q| 
1 

(2tt) 2 |Q| 



Jo 



dk kf-F{k) 
rfifc k 2 F(k) 



I 

Jo 

I" 

Jo 



1 



^0 

Q 2 



/c c (fc, w ) 



vfc 2 



2 \ 



1 + -4 

Q 2 



K c (k, w ) k> 



t/wo X 

W() 



(16) 



vk 2 



Q 2 v 2 /^ + a> 2 



cr c (k, w ) 



(2tt) 2 |Q| 
1 

(2n) 2 \Cl\ 



dk k z F(k) 



dk k z F(k) 



X" 

Jo 

I" 

Jo 



afwo 



afwo 



Q 2 



,2\ 



1 + 



Q 2 



X(k, wo) 



I(k, w ) 



and vanishing transport of angular momentum: (u x u y ) = 
(u x u y ) = 0. We now examine the dependence of the turbulence 
amplitude and transport of particles in Eq. ( fT~6b on the correla- 
tion time of the forcing C(r). 
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4.1. Infinitely correlated turbulence 

When C(t) = 1, the correlation functions in Eqs. (O and ( fTTT i 
take the following forms: 



1 



and cr c = . (17) 



^ Dk 2 ' Kc v 2 k 4 + aj 2 Q a " U " c vW + ^l 

Using Eq. (fT7l > in Eq. ( TToT l and performing the integration on 
the (l>o variable, we obtain the following equations: 



F(k) 1 



2 f + °° F(k) II O 2 + 1 i 

/ 1 2Q 2 - 5 

2 H r- h arctan Q, 

\ 2(1 +Q 2 ) 2Q„ 

' +0 ° F(Jt) 1 / 0» + l 

t/A: — — ~ 1-1 H — arctan Q, 



(18) 



XX 2 f 

n yy _ J_ f + °° J( 1 
r (2?r) 2 Jo 



Dv 2 k* £ll \ Q, 
'•' '• • ' ' il- - 1 



£>v 2 Ar 4 C2 2 



1 + 



a 



arctan Q t , 



where Q« = £2. (A) = \£l\/(vk 2 ). In the limit Q» 



the 



turbulent transport and intensity, both for the parallel and per- 
pend icular componen t s, all t end to zero as Q 1 . These agree 
with Kic hatinov et al.l (|1994). The effect of rotation on turbu- 
lence amplitude appears with Q», which becomes large on large 
scales. In other words, turbulence quenching is more severe on 
large scales, leading to effectively strong turbulence on small 
scales. This is consistent with the reduction of scale of motion 



as evidenced in numerical simulations (iBrummell et al. 1998). 

Interestingly, the turbulent transport of particles is propor- 
tional to D l , suggesting that the turbulent transport of heat 
scales as k with the molecular heat diffusivity k. Thus, in the 
Sun where k » D, the transport of particles is expected to 
be faster than that of heat by a factor k/D ~ 10 7 /10 2 ~ 10 s . 
However, this holds only for an incompressible fluid, as the 
equation for the transport of he at is not the same as that o f par- 
ticles for a compressible fluid dSpiegel & Veronis 1960). The 
relevance of this result for the solar context may thus be ques- 
tionnable. 



4.2. 6-correlated turbulence 

We now consider a short correlated turbulence modelled by 
C(t) = t/(5(t). In this case, the temporal correlations in Eqs. 
(0 and (fTTb become: 



LL 

2 



and 



. 



(19) 



leading to the following turbulent intensity and transport: 

2 r°° F(k)T f 



(4) = {u]) 



l r 



2 r 

2^ Jo 



WW Jo v 
F(k)T f 1 
v 2 k 2 Q 2 



dk 



(20) 



nl + d+b) 2 (a. 

-(1 + b) + arctan — 

2fi« \ 2 



dk . 



m = 



i r° F( 
i^Jo ~ 



F(k)Tf 1 



2(2. 
1 +b + 



o 



v 2 k 2 ill 
d+b) 2 



2Q, 



arctan I — 



/O, 
2 



dk 



Here, b = 1/Pr. Eq. (F20l shows that, in the case of a short 
correlated forcing, only the turbulent transport is suppressed 
(by a factor £2~') whereas the turbulent intensity is of the same 
order as that without rotation. This is because, in this case, in- 
ertial waves only have an effect on the phase of the velocity 
field and thus do not modify its amplitude. Alternatively, this 
is because the flow driven by a forcing with a short correlation 
time Tf < Q. 1 has no coherent motion to be affected by inertial 
waves (or rotation). This is to be compared with the result ob- 
tained in the case of an infinitely correlated forcing (Sect. 14. Il l 
where the turbulence amplitude was reduced as Q~' due to ro- 
tation. Furthermore, the dependence of the turbulent diffusivity 
on the molecular diffusivity is not as simple as in the infinitely 
correlated case. However, in the large rotation limit (Q* » 1), 
we can easily see from Eq. ( 1201 that the turbulent diffusivity 
depends neither on P r nor on D. Consequently, the transport of 
heat and particle is the same for rapid rotation. 

4.3. Finite correlation time with exponential correlation 
function 



In Sect. 14. II and 14.21 we considered two extreme limits where 
the correlation time was infinite or zero, respectively. Here, we 
discuss a more realistic case where the turbulence has a finite 
correlation time r y by using an exponential correlation function 
C(t) = exp[-2r/r/]. In this case, the functions in Eqs. $6fy and 
(fTTb can be simplified as: 



K c 



vk 2 + 2/Tf 



tc = 



(vk 2 + 2/T/) 2 +w 2 

OJO 

{vk 2 + 2/t/) 2 + u>l 
1 



(21) 



Dk 2 +2/r f ' 

Putting Eq. ( f2Tb in Eq. ( fToT ). we can calculate the turbulent in- 
tensity and transport as previously. We omit, however, the gen- 
eral form of the results which are too complex, and only pro- 
vide a result that is valid in the limit Q. — > oo. To this end, we 
can simply take the limit Q. — > oo in Eq. ( TToT l and obtain the 
following result: 



<"a> ~ ,„ 2 9,^, f dkk 2 F(k) [ da) Kc ^° ) , (22) 



duo 



dkk 2 F(k) I 
o Jo 

< M > 2) ~ TT-oH dkk 2 F(k) 

(2tt) 2 |Q| Jo Jo 

D t ~ TT-alFt dkk 2 F(k) du I(k, 
(27!\r|i2| Jo Jo 

°? ~ T^TToiFTi dkk 2 F(k) doj I(k, 
(2n) 2 \Q,\ Jo Jo 

We note that, in order to obtain Eq. ( l22l i. the limit Q. — > oo 
was taken inside the integral as the resulting integrals are all 



vk 2 ' 

*c(k, up) 
vk 2 ' 

(jJq) . 
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convergent. Eq. ( l22b shows that both turbulent amplitude and 
transport scale as Qr l for large rotation rate. This is similar to 
the case of the infinite correlation time case. Furthermore, Eq. 
( l22l shows that the ratio between the parallel and perpendicular 
transport of particles (or heat) is given by: 



<"?> 



D X T 



= 2 for a 



(23) 



again similar to the case of C(r) = 1. We note that the large 
rotation limit Q — > oo taken above requires that not only Q. » 
vk 2 but also Q » tZ 1 . In other words, the correlation time 
of the forcing should be larger than Q -1 , In the solar context, 
Q ~ 3 x 10 _6 i. Thus Tf must be larger than a few days. When 
Q < Tj l , it can easily be shown that the results become similar 
to those in the case of ^-correlated forcing studied in Sect. 14.21 
(i.e. (u 2 ) ~ Q° and D T ~ Q" 1 ). 



4.4. Power law correlation 

In view of the highly intermittent nature of stellar activity, it 
is interesting to consider the case with a correlation function 
of the forcing as a power-law. For simplicity, we consider a 
correlation function given by: 



(24) 



with < n < 1 (note that it is also possible to consider a func- 
tion with exponent fi > 1 by introducing a cutoff frequency). 
To calculate the various correlation functions, we need the fol- 
lowing integral: 



Jo 



(vk 2 ) I 



dt , (25) 



where the last equation is obta ined using steepest de scent 
method. Using Watson's lemma (Ben der & Orszagll 19751) . the 
asymptotic behaviour of this integral for (Dq = ojo/vk 2 » 1 can 
be found as: 



m-n) — 

Wo 



(26) 



Thus the correlation functions in Eqs. © and (fTTT i become: 

_ rq-M) 

ic ~ Dk 2 f*' 



(27) 



(Tc 



vk 2 
vk 2 



(l-j"); 



where t/» = Tjvk 2 and «o* = oJo/vk 2 are the correlation time 
and the oscillation frequency of the noise scaled by the diffu- 
sion time, respectively. Putting Eq. (|27T i in Eq. ( TToT i. we then ob- 
tain the turbulence amplitude and turbulent transport of chem- 
icals for an isotropic forcing with a power-law correlation: 



(2 



2 r + 

W Jo 



dk 



F(kV ft 
v 2 k 2 Q t 



cos 



(1 



2Qt 



2\u(u + 2) 



(28) 



Dt = 



1 



(2tt) 2 
2 



f 

Jo 



(2n) 2 I 

' X* 

Jo 



dk 



dk 



dk 



— C' 

y 2 £ 2 Q„ 

Dv 2 /t 4 Q» ' 
F(k)Tf, , 



0-4 



2Qi + ljOg 
HiH + 2) 



(2n) 2 J,, Dv 2 k 4 D. t ' 

in the rapid rotation limit f2» = |Q|/(vfc 2 ) » 1. T is a con- 
stant which is independent of Q. From Eq. (f28]l, we can see that 
the turbulence intensity scales as l^ -1 for rapid rotation. This 
shows that the scaling of the turbulent intensity depends sen- 
sitively on the exponent of the power law. Note that the result 
of infinitely correlated turbulence is recovered for fi = 0. In 
comparison, the turbulent diffusivity in both directions scales 
as |Q|~' independent of the power exponent fi. Furthermore, 
we can see that the ratio between the turbulence in the x- and 
y-direction depends only on the power exponent: 

<^> 
<«?> 



jU + 1 



(29) 



Consequently, the turbulence is always more vigorous in the x- 
direction (along the rotation axis) than in the y-direction. We 
recover the factor 2 for /i = 0, corresponding to the exponen- 
tially correlated case, see Eq. (1231 . As the power exponent is 
increased, the turbulence becomes more isotropic. Note that the 
ratio between the transport of species in the x- and y-direction 
is always 2, as in the case of an exponential correlation func- 
tion. 

To summarise, for an isotropic forcing and rapid rotation, 
we showed that the turbulent diffusivity is reduced as QT l , 
irrespective of the correlation time of the forcing, with the 
transport in the direction parallel to the rotation vector be- 
ing twice as fast as than in the direction perpendicular to 
the rotation. This would lead to a warmer pole than equator 
dKichatinov & Riidigerll 19951) . In comparison, the turbulent in- 
tensity is found to be reduced with different scaling with O 
depending on the properties of the correlation function of the 
forcing. Specifically, in the case of a finite correlation time, the 
turbulence intensity scales as CI 1 for rapid rotation, while it 
is independent of the rotation rate for a short-correlated forc- 
ing (Sect. I4.21 i. Furthermore, for a temporal correlation func- 
tion given by a power-law, turbulence amplitude scales as Q." 
with -1 < n < 0. 

5. Anisotropic forcing 

We now consider the case where the forcing f is anisotropic 
as should be relevant for the turbulence in the convection zone 
(for instance, due to the underlying stratification). We are par- 
ticularly interested in the Reynolds stress as the combination 
of anisotropy and rotation can lead to non-diffusive fluxes of 
angular momentum. In the case where g is a unit vector in the 
direction of an isotropy, the c orrelation function of the forcing 
can be written ( Riidigerll 1989b : 

</Kk, f)/)(k', t')) = (2;r) 3 C(\t - I'D 6(k + k') G(k) x (30) 

(g-k) 2 g-k 
— 2 —Sij ~ 8i8j + —(8<kj + gjkd . 
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This forcing in Eq. ( f30l > is highly anisotropic with only a non- 
vanishing component perpendicular to g: 



<(f ■ g) 2 ) = , 



(31) 



<(fxg) 2 ) 



1 r 

' 2^ Jo 



G(k)k 2 dk . 



In the following, we assume g to be in the vertical direction. 
In general, one can consider a less anisotropic forcing by com- 
bining the isotopic part of the forcing (TTZl l and anisotropic part 

CP. 

5.1. Rotation parallel to the vertical direction 

When the rotation vector is parallel to the vertical direction i.e. 
g || Q (as in the case of the turbulence near the pole), Eq. (l30b 
translates to: 



= 012 = and <p22 = 



r 



y + a 



jG(k) 



(32) 



Using Eq. (l32l in Eqs. © and (O, we can show that the turbu- 
lence amplitude becomes 



(2/r) 2 |Q| 
1 



dkk 2 G(k) 
o Jo 



,\Q\ 

dkk l G(k) I dtoo 
Jo 

vk 2 ^ 



2\ 



2vk 2 2(v 2 k 4 + cj 2 ) ) 
1 

(2n) 2 \Q\ J(1 



i _ ±5* 
n 2 



0J() 



2(v 2 k 4 + wg) 



/-co wn| 
dkk 2 G(k) duo 
Jo Jo 



2\ 



n 2 



(33) 



1 + 



,,2\ 



Q 2 



vk 2 K c — u>oO~ c 



2vk 2 2(v 2 k A + w 2 ) 



while the transport of angular momentum vanishes. Here again, 
we can see by using Eq. d2TT > that all the integrals in Eq. ( 1331 
converge in the limit Q — > oo. Thus, turbulence amplitudes are 
reduced as Q -1 for large rotation rate. We do not provide here 
the results for the turbulent diffusivities as they are very similar 
to those derived in the isotopic case. Thus, in the rapid rotation 
limit, turbulent diffusivities scale as Q 1 for all the different 
types of correlation function of the forcing that were considered 
previously. 

5.2. Rotation perpendicular to the vertical direction 

Near the equator, the rotation vector is perpendicular to the 
vertical direction. Since the rotation is assumed to be in the 
x-direction, the vertical direction for anisotropy can be chosen 
to be in the z-direction without loss of generality. Eq. d30i > then 
translates to: 



0ii =(y + a 2 )G(k) 
and < 



'22 



»ia = a/3G(k) 
a 2 /3 2 



(34) 



y + a L 



T G(k) . 



In this case, the turbulent amplitude and diffusivities can 
be shown to be very similar to those derived in the isotopic 
case. Consequently, the scalings are the same as those derived 



in Sect. [4] However, in this case, the transport of angular mo- 
mentum does not vanish as discussed below. Putting Eq. ( 1341 
into Eq. ([8]) and keeping only the even terms in ft and a (as 
the odd terms vanish after angular integration), straightforward 
but cumbersome calculations lead to the following expressions 
for the correlation function of the y and z components of the 
velocity: 



m v m,) = \ d^k 

2(2tt) 3 J 



9a(l + a 2 )G(k) o>qk c + vk 2 cr c 



(y + a 2 ) 3 / 2 



v 2 k 4 + a> 2 



(35) 



while the other correlation functions vanish. Performing the in- 
tegration over the angular variables, Eq. d35l l can be simplified 

as: 



x C°° f |n| 

(u y u z ) = - dkk G(k) t 

2(2nYLl l Jo Jo 



da>o x 



(36) 



wo 



1 + 



Q 2 



CJqKc + Vk (T c 



v 2 k A - 



where x 1S tne sig n °f The Reynolds stress in Eq. ([35} 
depends on the direction of the rotation, re presenting non - 
diffusive momentum transport, i.e., A effect dRudigerlll989l) . 
This effect favours the creation of velocity gradient rather than 
smoothing it out and can thus be responsible for the differential 
rotation of stars. 

We now examine the behaviour of Eq. ( |35t for different 
temporal correlation functions. In the case of an infinitely cor- 
related forcing, Eq. ( |36l i can be shown to become: 



(2^1 



2(2. 
3 +H 2 



dk 



G(k) 

v 2 k 2 nl 



(37) 



^ 2 -3 
H — arctan £1, 



1 + O 2 Q. 

Eq. ( T37l l shows that the non-diffusive fluxes of angular momen- 
tum scale as O" 2 in the large rotation limit. This agrees with 
iKichatinovl (1198a) . Similarly, we can show that for a finite cor- 
relation time (with exponential form), the limit Q.* — > oo can 
be taken inside the integral of Eq. d36b and consequently, the 
non-diffusive fluxes of angular momentum scale also as Q" 2 in 
the large rotation limit. 

In comparison, when the forcing is (5-correlated, Eq. 
becomes: 



TfG(k) 



u v u 7 ) — I dk — — x 

> J 2(2n) 2 Jo 6vQ4 

[Q„(4Q 2 - 3) + 3(1 - Q 2 ) arctan n,] 



(38) 



which shows that the non-diffusive fluxes of angular momen- 
tum scale as Q~' in the rapid rotation limit. 

Finally, for a power-law correlation function, by putting Eq. 
d27l > into Eq. d36l l. we obtain the non-diffusive fluxes of angular 
momentum as: 



(u y u z ) 



W Jo 



dk 1 - - r(i 



2(27t) 2 J """ vk 2 nl 
OH cos [(l-n)n/2] 
2(// + 2)cos [(1 +h)tt/2] 



cos 



■jU)X 



(l + 4 



(39) 
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in the rapid rotation limit. In this limit, the A-effect scales as 
D.^ 2 . Thus, for a power-law correlation function, the A-effect 
scales as Q." with 1 < n < 2. 

In summary, we show that non-diffusive fluxes of angular 
momentum do not vanish only when the rotation is perpendic- 
ular to the inhomogeneity (such as stratification in the radial 
direction or shear in the azimuthal direction). This is in agree- 
ment with previous results. Furthermore, we show that these 
non-diffusive fluxes scale as Q 2 - Q -1 , in the rapid rotation 
limit, depending on the property of the temporal correlation 
function of the forcing. 

6. Multi-scale renormalisation group analysis 

The results presented in the previous sections are derived us- 
ing the quasi-linear theory. This is strictly valid only for two- 
scale turbulence, with a spatial gap between a large-scale for 
the mean-flow and a small scale for the fluctuations, and for 
weak turbulence. Therefore, the applicability of this approxi- 
mation is limited. In this section, w e utilise a renormalisation 
group technique (see iMoffattll 1 98 ll and references therein for 



details) to consider a more realistic situation of a strong turbu- 
lence with a no-clear scale separation between turbulence and 
mean fields. This method is based on the computation of the 
turbulent diffusivity iteratively by considering a range of scales 
lis with: 



l n > h-\ > ■ ■ ■ > h > 



(40) 



and /o is the scale at which the molecular diffusion becomes 
dominant. The effective turbulent diffusivity D n+ \ for scale n+ 1 
is then calculated by using the value of the turbulent diffusiv- 
ity at the preceding step, i.e., D„ via quasi-linear theory. Since 
the turbulent diffusivity scales as D l [see Eq. (1181) 1. we can 
calculate D\ in terms of Do as: 



1 pVk 

Do Jm 



P(k, v, Q) dk , 



(41) 



where P(k, v, Q.) is a function depending only on modulus of 
the wave number, molecular viscosity, and rotation. Note that 
the integration range is over the scale I such that Iq < I < l\. 
Repeating this procedure iteratively, one can compute the tur- 
bulent diffusivity at scale n + 1 as follows: 



i r l »» 
= -r — ^ P(k,v,n)dk. 

Do + D„ J 1//n+1 



(42) 



If we take the limit n — > oo, we obtain a continuous range of 
scale and turbulent diffusivity. Thus, writing D(k) as the turbu- 
lent diffusivity acting on scales smaller than k , Eq. d42l) can 
be rewritten 



dD(k) = 



1 



Do + D(k) 



P(k) dk 



(43) 



The turbulent diffusivity Dj is obtained by integrating over all 
scales in the inertial range, starting from n — with the value 
of molecular diffusivity Do = D, with the result: 



Di 



P(k, v, Q) dk 



1/2 



(44) 



Eq. ( PPfl i shows that if the quasi-linear theory gives a turbu- 
lent diffusivity as Q. ", then the renormalisation analysis pre- 
dicts a scaling Q.~ n ^ 2 . Consequently, as the quasi-linear theory 
predicted a scaling as Q -1 (irrespective of the form of tempo- 
ral correlations), the turbulent transport of particles (or heat) 
scales as Q~'^ 2 in the large rotation limit. Interestingly, this re- 
sult has a weaker dependence on Q compared to that according 
to quasi-linear theory. Furthermore Z>" /Dj — V2 in the rapid 
rotation limit, with a weaker anisotropy compared to the quasi- 
linear result [see Eq. (l23l)l. 

Finally, we note that quasi-linear or multi-scale renor- 
malisation group analysis cannot capture a strong anisotropy 
in the rapid rotati on limit with a ten dency towards two- 
dimensionalisation dCambon et al.l 19971) . In order to capture 
the highly anisotropic nature of rotating turbulence, one needs 
to use more sophisticated closure, such as wav e- turbulence 
which leads to anisotropic spectra dGaltierl 12003). In the ex- 
treme limit of rapid rotation, this rotation-induced anisotropy 
could be captu red by assumin g geostrophic balance in the equi- 
librium state (Pedlovsky 1983). 



7. Summary and discussion 

We studied the structure of turbulence in a rotating medium 
with an arbitrary external forcing. Using the quasi-linear the- 
ory, we first performed a thorough study of the turbulence in- 
tensity and transport (of chemical species or heat) in the case 
of an isotropic forcing by considering different types of tem- 
poral correlation functions. We also examined the transport of 
angular momentum in the case of an anisotropic forcing. 

Specifically, we showed that for an isotropic forcing, the 
transport parallel and perpendicular to the rotation vector has 
the same scaling with the rotation frequency £2: they are both 
reduced by a factor Q 1 compared to the case without rotation 
irrespective of the temporal correlation function of the forcing. 
Furthermore, the transport in the direction parallel to the rota- 
tion is twice as fast than the one in the perpendicular direction 
[see Eq. j23l l. This result was shown to be robust independent 
of the specific form of temporal correlation C(t) (exponential 
and power law) or correlation time (finite , short or infini te). 
Note that this result was shown by iKichatinov et alJ (119941) in 
the case of an infinitely correlated forcing. 

In comparison, the scaling of the turbulence intensity de- 
pends on the property of the forcing. For exponential corre- 
lation function with infinite memory, all the components of 
the velocity amplitude scale as Q~' while their intensity in 
the direction parallel to the rotation is twice that in the per- 
pendicular direction (similar to the transport). However, for 5- 
correlated turbulence, the turbulence amplitude becomes inde- 
pendent of the rotation rate. Thus, in this case, inertial waves 
affect only the cross-phase of the velocity field and not its am- 
plitude. Furthermore, in the case of a finite correlation time 
with a power-law, the turbulence intensity scales as Q" with 
-1 < n < while the ratio of the turbulence amplitude in the 
direction parallel to that perpendicular to the rotation can be 
any number between 1 and 2 depending on the exponent of the 
power-law. 
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In the case where the driving force is highly anisotropic, 
it was shown to give rise to non-diffusive fluxes when the 
rotation is perpendicular to the anisotropy. In the large rota- 
tion limit, the transport of angular momentum scales as CI' 2 
for a forcing with a finite correlation time. However, for a 6- 
correlated turbulence, we found that the (quenching) effect of 
rotation was reduced with non-diffusive fluxes now scaling as 
Q . Furthermore, for a power-law correlation function, the 
momentum flux was found to scale as Cl p with — 2 < p < — 1. 

Finally, we performed a multi-scale renormalisation anal- 
ysis which permitted us to go beyond the quasi-linear theory 
in the calculation of the turbulent diffusivity. We found that, 
in comparison with the quasi-linear theory result of a turbu- 
lent diffusivity proportional to Cl~ l , the renormalisation theory 
predicts a scaling Q ^ 2 . Consequently, we expect the turbulent 
diffusivity to depend weakly on CI. It also gave us an anisotropy 
in transport weaker than the quasi-linear prediction, the parallel 
component being V2 faster than that of the perpendicular. 

We should however note that the scaling of our results can 
be modified in the case of a bounded domain of finite size. 
In the calculation of all the turbulent transport coefficients, we 
obtained a result proportional to the following type of integral: 



/(k, Q) 



-/ 



H(k) 
v 2 * 4 + col 



(45) 



where coo - {CI- k)/k is the projection of the unit vector in the 
direction of the wave number on the rotation axis. When the do- 
main of integration is unbounded, the integration over the an- 
gular variable makes this integral proportional to Q 1 , when the 
rotation rate CI is sufficiently large [e.g. Eq. (I18ll l. This is be- 
cause this integral involves some contribution of order 1 (when 
CI ■ k = 0) and others of order Q~ 2 . However, in a realistic situ- 
ation, the domain of integration (in Fourier space) is bounded. 
That is, there is a minimal wavenumber (corresponding to a 
maximum length, for instance the size of the box) in the direc- 
tion of the rotation, for example, k m = min(fe A ). The aforemen- 
tioned scaling of Q _1 is valid only when v 2 k 6 » Cl 2 kf n . In the 
opposite case (v 2 k 6 «: Cl?k 2 n ), the term coi in Eq. (05]) is always 
dominant, giving the scaling of CI 2 for large rotation rate. In 
this case, the multi-scale analysis would give a scaling of Cl l 
instead of Q" 1 ^ 2 . 

8. Implications for stellar convection zone 

We discuss here some of the important implications of our re- 
sults for the Sun and other stars. For parameters typical of the 
Sun, CI ~ 2.3 x KTV 1 and v ~ 10 2 cmV\ the effect of 
a finite-size domain discussed in Sect. [7] can be ignored if the 
characteristic length scale is larger than (v/Q) 1 ^ 2 ~ 1.7x1 4 cm. 
This is likely to be the case in the solar convection zone as 
convective motions can extend to a distance of approximately 
4 X 10 9 cm. In this case, the quasi-linear theory predicts that the 
turbulent diffusivity scales as Q _1 for rapid rotation while the 
turbulence amplitude scales as Cl°-Cl l depending on the prop- 
erty of the forcing. In comparison, the renormalisation analysis 
predicts a turbulent diffusivity scaling as Cl~ l l 2 . This result thus 
suggests that the effect of rotation o n particle or heat trans port 
is weaker than previously thought dKichatinov et alJll994ft . If 



the effect of rotation had been much stronger in reducing trans- 
port, it could have inhibited the mixing of light elements in the 
convection zone, th ereby leading to a weak mixing as required 
dBarnes et alj [l999). Our results have confirmed that this can- 
not possibly be true and that rotation cannot play an important 
role in reducing the mixing of light elements in the solar con- 
vection zone to explain the o bserved (rather modest) dep letion 
of light elements on the Sun (Schatzman & Baglin ll99ll) . 

Traditional stellar modelling does not involve any 
anisotropy in the turbulent diffusivity of particles (or equiva- 
lently in the turbulent conductivity of temperature) due to ro- 
tation. We have shown here that such an anisotropy can be in- 
duced by rotation. This is however not a significant effect as 
the transport in the direction parallel to the rotation is only 
twice as fast as that in the perpendicular direction, according to 
the quasi-linear theory. This anisotropy is quite robust against 
the temporal correlations used and ha s been shown to give 
rise to a warmer pole than the equator dKichatinov & Riidigei 
Il995h . This anisotropy is however shown to be weaker ac- 
cording to the multi-scale renormalisation analysis, the ra- 
tio being only V2. Consequently, rotation may not be suffi- 
cient for inducing anisotropic transport in stars, for instance, 
in order to cause a sufficient latitudinal temperature gradient 
in the convective envelopes as required by numerica l models 
dKichatinov & Rudiger|[l995l: iKiiker & Rudigerll2005l) . As we 
have shown in our previous studies, an anisotropic transport 
can easily be in duced by a strong shear layer (such as the solar 
tachoc line, see Kim 2005; Leprovost & Kim 2006), magnetic 
fields dKiml l2006t iLeprovost & Kimll2007at [Kim & Leprovostl 



2007ah or stratification (Kim & Leprovostl Fo07bl) . In these 
works, shear flows are also shown to quench turbulent trans- 
port significantly leading to weak mixing. 

Finally, we note that our results depend on the magnitude 
of Q, the diffusion rate vk 2 (or Dk 2 ) and the characteristic fre- 
quency of the forcing 1/t/. Obviously, stars of different spec- 
tral types have different values for CI, v, D and tj. It is thus 
important to understand the origin of tj and to perform a sys- 
tematic study on the development of turbulence theory as a 
function of masses and ages. This will be addressed in future 
studies. 
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Appendix A: Turbulence amplitude and transport 

In this appendix, we provide the form of A'j k (k, too) and 



M' jk (k, coo) used in Eqs. © and ([8]t. 
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Here again, the functions k c , ct c and ( c characterising the in- 
fluence of the temporal correlation of the forcing on the turbu- 
lence are defined by Eqs. © and (fTTT >. 
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